Average shape of fluctuations for subdiffusive walks 
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We study the average shape of fluctuations for subdiffusive processes, i.e., processes with uncor- 
rected increments but where the waiting time distribution has a broad power-law tail. This shape 
is obtained analytically by means of a fractional diffusion approach. We find that, in contrast with 
processes where the waiting time between increments has finite variance, the fluctuation shape is no 
longer a semicircle: it tends to adopt a table-like form as the subdiffusive character of the process 
increases. The theoretical predictions are compared with numerical simulation results. 

PACS numbers: 05.40.-a, 05.45.Tp, 45.10.Hj 
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I. INTRODUCTION 

Complex systems are often described by their lack of 
a characteristic length or time scale over many orders 
of magnitude, which gives rise to events whose distri- 
bution in sizes is a power law with no characteristic 
size (fractal behaviour). Examples are everywhere 0: 
seismic activity, turbulence, solar flares, Brownian mo- 
tion, length of rivers and blood vessels,. . . In particular, 
a power-law distribution with no characteristic tempo- 
ral size events appears in the analysis of stock price 
changes 0, Q- 

river floods g , Barkhausen noise jE@, 
glassy systems 0, H, 0, atomic cooling 0, and fluo- 
rescence of quantum dots 0, I* 1 these cases the 
resulting dynamics is strongly intermittent, with bursts 
of activity separated by long quiescent intervals. 

When these temporal intervals are waiting times St 
between jumps of size Sx, then the stochastic process x(t) 
can be seen as the trajectory of a subdiffusive random 
walker (provided that the variance of 5x is finite). A 
typical subdiffusive trajectory is shown in Fig.^ Systems 
that exhibit anomalous subdiffusion characterized by an 
anomalous Fick's second law 



2K-, 



r(i + 7) 



(i) 



where < 7 < 1, are ubiquitous in nature 0,0]. [K 1 is 
the (generalized) diffusion constant and 7 is the anoma- 
lous diffusion exponent.] Also, models based on subdif- 
fusive random walkers are useful for understanding com- 
plex systems. Two nice examples are the "trap model" , 
proposed to explain aging in disordered systems H, , 
and the comb model, to understand diffusion phenomena 
in complex structures such as percolation clusters [15j . 

Recently H B E E3 the study of the average shape 
of the fluctuations of stochastic processes x{t) has been 
considered as a useful tool to gain insight into the sys- 
tem that generates the fluctuation. Thus, it has been 
argued |18| that the average shape of fluctuations is a 



better tool for discriminating between theories than crit- 
ical exponents. 

In Ref. [lj| Baldassarri et al. consider the average 
shape for a stochastic process of the form x(t+l) = x(t) + 
5x, where Sx is a random variable. Let us denote the 
average shape of fluctuations of time span T by (x(t))T- 
Baldassarri et al. find that the shape follows the scaling 
law 



(x(t)) T =T<*/ 2 f(t/T), 



(2) 



where a is the diffusion exponent and / is a semicircle 
whenever 5x follows a distribution with finite variance 
(Gaussian walks) or a distribution X(Sx) with a broad 
power-law tail, X(Sx) ~ (Sx)^^ 1 with < fi < 2, so 
that the variance is infinite (Levy flights). This result had 
already been obtained by Fisher for Gaussian walks 0, 
Sec. 7.1]. The fact that the average shape of fluctuations 
is a semicircle for both Gaussian walks and Levy flights is 
a nicely surprising result that led us to wonder to what 
extent it might hold for other walks with uncorrelated 
jumps. In particular, we investigated the average shape 
of subdiffusive stochastic processes (subdiffusive walks) 
of the form 



x(t + St) = x(t) + Sx 



(3) 



where Sx are uncorrelated increments that follow a dis- 
tribution with finite variance and St is a random vari- 
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FIG. 1: Subdiffusive trajectory with Sx = ±1 and 7 = 0.9. 
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able whose distribution tp(St) has a broad power-law tail: 
ip(5t) ~ with < 7 < 1. Figure [J shows the 

position x{t) for a stochastic process of this kind with 
7 = 0.9 and with jumps of unit length (8 = ±1) made 
with equal probabilities. 



II. AVERAGE FORM OF A SUBDIFFUSIVE 
FLUCTUATION 



The average shape of a fluctuation can be calculated 
through the expression 0, 0] 

/°°dxxft(x,i|{x(0)};xo,0,x o ,:r) 

w))t = lim " — (4) 

x ->o+ J efx S2(x, r|{x(0)j; xq, 0, xq, 1 ) 



as [13|: 

P(x,t) 



ttW 

H ll 



(1-7/2,7/2) 
(0,1) 



(6) 

Taking into account that the Laplace transform of P{x, t) 
is 

P(x,u) = ^=exp(-Ju'r/K J \x\) (7) 



one finds for x > 

u 7/2-l 



F(x,u|x ,0) 



{ e -° (a:o - K) [e(x)-e(x-xo)] 

o(x+x )\ ^ 



+e- a( ^ a;o) e(x-xo 



which implies 
F(x,t\x -> 0,0) 



where 



where f2(x, t|{x(0)}; Xo, 0, xo, T) is the probability that 
the walker with trajectory {x(0)} for t < is at x at 
time i provided that he was at xo > at time and at 
xo at time T without ever touching the axis at x = 
in the time interval (0,T). The upper line in Eq. Q 
means average over all the trajectories {x(0)} that reach 
x at time t = 0. Let F(x,t\{x(t')};x',t') be the prob- 
ability that the walker with trajectory {x(t')} for t < t' 
and who was at x' > at time t' reaches x > at 
time t without ever touching the axis at x = 0; let 
f2(x, i|xo, 0, xo, T) be the probability that the walker is 
at x at time t provided that he was at xo > at time 
and at xo at time T without ever touching the axis at 
x = in the time interval (0,T); and let F(x, t\x',t') be 
the probability that the walker who was at x' at time 
t' reaches x > at time t without ever touching the 
axis at x = 0. For walks without memory (Markovian 
walks) J7(x, t|{x(0)}; xo, 0, xo, T) = f2(x, t\xo, 0, xo, T), 
F(x, t\{x(t')}; x', t') — F(x, t\x', t'), and one can write 
f2(x, t|x , 0, xo, T)as Fix, t\xo, Q)F(x, T—t\xo, 0), so that 
Eq. J2J becomes p^.ll7| 

dxx F(x,t\x o ,0)F(x,T - t\x ,0) 

{x[t))T = hm . 

-o-o+ f™dxF(x,t\x ,0)F(x,T-t\x ,0) gM/T) = 

(5) 

This equation is not exact for subdiffusve walks because 
they are not Markovian. However, for subdiffusive walks, 
the memory, i.e., the effect of the fact that at t' the parti- 
cle was at x' on the probability that the particle at time 
t > t' is at x, decays as (t - i')-T/r(l - 7) 13}. This 
implies that the approximation of F(x, t\{x(t')}; x', t') by 
F(x, t\x', t'), and, consequently, the accuracy of Eq. (J5J 
for subdiffusive walks, improves when t — t' increases and 
7 is close to 1. 

The probability F(x, t|xo,Q) can be calculated by 
means of the method of images F(x,t\xo,0) = P(x — 
Xo,t) — P(—x — Xo,t) [13,01, P(x — xo, t) being the prob- 
ability density that the free process (without boundary 
conditions) that at time t < was at x is at x at time 
t. For subdiffusive processes, and for t _7 /r(l — 7) <C 1 
|2^ , P(x, t) can be written in terms of Fox's H function 



where a = \Ju* j K 1 and O(x) is the Heaviside step func- 
tion. As we are interested in the limit xo — > with 
x > x , we expand the term inside the bracket in powers 
of xq and get 



F(x, u\xq — > 0, 0) = Xq 



7— 1 

U ^ c -Ju-</K-,x ^ 



X H 10 



(1-7,7/2) 
(0,1) 



(10) 

Inserting this expression into Eq. JSJ and carrying out 
the integrations |2^| one finds that the average shape 
of a subdiffusive stochastic process (subdiffusive random 
walk) is given by 



(x(t)) T - T^ 2 f 7 (t/T) - yfK^g^t/T) (11) 



I t TTll 

\TJ 22 


( t/T y/ 2 

\l-t/T) 


(-l,l),(l-7,i)" 

(0,1), (o,|) 


H 22 


( t/T \ 


7 /2 


(0,1), (I-7,?) " 

(0,1), (i,i) 





(12) 

In Fig. we plot the (normalized) average shape of fluc- 
tuations for several classes of subdiffusive processes. We 
see that the shape tends to a table-like form as 7 de- 
creases, i.e., as the subdiffusive character of the process 
increases. 

Of course, for 7 — > 1 one recovers the Gaussian result 
h{t/T) = y/l6D/iry/t/T(l - t/T) HHH because the 
upper and lower Fox's H functions in Ea. H12(l become 
2 z/tt (1 + z 2 ) 2 and z/2 y 7 ^ (1 + z 2 ) 3 / 2 , respectively. 

The area N(T) of the average fluctuation of duration 
T is given by 



N(T) = / ds(x{sT)} T 
Jo 



(13) 
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FIG. 2: Normalized average fluctuation for subdiffusive pro- 
cesses for several values of 7. The shape is normalized so that 
its area is 1, i.e., we plot (x{t))T /N{T). The lines are the 
theoretical result for (at t/T =1/2 and from top to bottom) 
7 = 1,0.9,0.75,0.5. Inset: simulation results for these same 
values of 7. 



TABLE I: The coefficient n 7 = dsg y (s) calculated by nu- 
merical integration. 



7 1/4 


1/3 


1/2 


2/3 


3/4 


4/5 


9/10 


1 


n 7 0.493 


0.501 


0.537 


0.612 


0.668 


0.708 


0.798 


v^/2 



where s = t/T. From Eq. Ullfl one finds 



N(T) = n^y/K^T-r 



(14) 



where n 7 = dsg^(s). Then, the normalized average 
fluctuation is given by 



(x(t)) T _ g 7 (t/T) 
N(T) n 7 



(15) 



We have not been able to calculate n 7 analytically. In 
Tabled we give some values evaluated numerically. 




FIG. 3: Normalized average fluctuation for the subdiffusive 
process with 7 = 0.9. The line is the theoretical result and the 
symbols are simulation results for T = 10 4 , 10 5 , 10 6 , 10 7 , 10 s . 




FIG. 4: Normalized average fluctuation for the subdiffu- 
sive process with 7 — 0.75. The line is the theoreti- 
cal result and the symbols are simulation results for T = 
10 4 ,10 5 ,10 6 ,10 7 ,10 8 . 




FIG. 5: Normalized average fluctuation for a process with 
7=1. The line is the semicircular form theoretical result and 
the symbols are simulation results for T = 10 3 (triangles), 10 5 
(squares), and 10 7 (circles). 



III. SIMULATION 

We carried out simulations of the fluctuations of the 
stochastic process Q where Sx takes the values +1 or 
— 1 with equal probabilities and where the waiting time 
5t between jumps follows the Pareto distribution ip(t) = 
7/(1 + <) 1+7 . In this case, the diffusion constant K 1 
is given by P3 K 1 = 1/[2T(1 - 7)]. The simulation 
results follow the pattern found in the precedent section: 
the fluctuation shape tends toward a table-like form as 7 
decreases (see inset in Fig. 

In Figs. |21 and Q] we compare the theoretical predictions 
and the simulation results for 7 = 0.9 and 7 = 0.75, re- 
spectively. The agreement is reasonable. We attribute 
the differences to, first, the approximate nature of Eq. (JSJ 
for non-Markovian processes and, second, to a slow con- 
vergence that requires longer T to set in. This can be 
clearly appreciated in Fig. where one sees that the well- 
established theoretical semicircular shape is approached, 
although very gradually, as T increases . It is even more 
gradual for smaller 7. We have not explored larger values 
for T because of the excessive computer time required. 
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FIG. 6: Area N(T) of the fluctuation for 7 = 0.75 (circles) and 
7 = 0.9 (squares). The lines are the corresponding theoretical 
prediction; see Eq. (1141 . 



with a power-law waiting time distribution ip(t) ~ i _1 ~ 7 
where < 7 < 1 (subdiffusive walks). Although the 
spatial increments Sx in this stochastic process are un- 
corrected, we find that, in contrast with Gaussian walks 
and Levy flights, the average shape of a fluctuation is no 
longer a semicircle: its form becomes flatter at the top 
and steeper at the extremes as the subdiffusive character 
of the process increases (i.e., as 7 decreases). 

Some possible sequels of the present work are obvious: 
one could consider the effect of correlated increments Sx 
(as done by Baldassarri et al. 16] ) or investigate the 
average shape of a fluctuation for Levy walks (in which 
\5x\ and St are proportional and follow a broad power-law 
tail). However, for these cases, the task of getting ana- 
lytical results, even approximate, such as those reported 
in the present paper certainly looks formidable. 



In the simulations we also calculated the area N(T) of 
the fluctuation. In Fig.|§]we plot N(T) for several values 
of T and for two values of 7. The agreement between 
simulation and theory is reasonable again. 

IV. CONCLUSIONS 

We have analyzed the average shape of the fluctuations 
of time series generated by the stochastic process Q 
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